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ABSTRACT 


The structure and behavior of barotropically unstable 
and stable wave disturbances in the vicinity of a zonally 
varying easterly jet are studied numerically with a 
linearized barotropic vorticity equation on a 8-plane. The 
easterly jet is approximated by a Bickley jet with a slow 
Zonal variation. The numerical results are also compared 
with a simple mechanistic analytical model using the local 
phase speed and growth rate concepts. The results are 
grossly similar in several respects to that expected from 
the parallel flow theory of barotropic instability, however, 
the resultant structure of the waves causes a spatial growth 
rate greater than predicted by the local growth rates 
computed with a parallel flow model. In the stable region, 
the structure leads to strong dynamic damping. When a 
uniform advective velocity is added to a variable mean flow, 
the differences between the behavior of the computed waves 
and that implied by the parallel flow theory are somewhat 
reduced. The waves remove kinetic energy from the mean 
flow and most of this energy is removed on the downwind side 
of the jet. The computed structure and behavior of the waves 
have a number of features that resemble those observed in 
the vicinity of the upper troposphere easterly jet during 


the summer monsoon. 
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1. INTRODUCTIQN 


In the Northern Hemisphere summer, a moderately strong 
easterly jet exists south of the Tibetan high near the 200 mb 
level (Krishnamurti, 197la, 1971b). Synoptic-scale moving 
disturbances occur at the level of the jet, and it appears 
that these disturbances arise from barotropic instability of 
the mean flow. The jet contains regions of large vorticity 
gradients where the necessary condition for barotropic 
instability is sometimes locally satisfied. If the observed 
disturbances were the result of barotropic instability, they 
would extract energy from the mean zonal flow and the plane- 
tary-scale waves, since the latter combine with the zonal 
flow to give the large vorticity gradients south of the 
Tibetan high. In fact, Kanamitsu et al. (1972) have shown 
that wavenumbers 6-8 in the wind spectrum in the region 
between 15S and 15N receive energy through barotropic 
interaction with the zonal and wavenumber ] flow. 
Furthermore, Krishnamurti (1971a) studied the wavenumber 
spectra of the meridional wind component for selected 
tropical latitudes and he found a peak in the spectra near 
wavenumbers 6-8 at latitudes near the easterly jet. 

Colton (1973) studied barotropic interactions between 
quasi-stationary long waves and transient synoptic waves 
using a semi-spectral numerical model. His long waves were 


forced with a specified divergence field following the 
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diagnostic model of Holton and Colton (1972). Results from 
a long term integration of his model reproduced several 
features of the upper tropospheric general circulation. In 
particular, his model simulated the disturbances observed by 
Krishnamurti (1971a, 1971b) in the vicinity of the easterly 
jet. His model showed that wave disturbances entering 
upstream of the easterly jet regime grow and increase their 
speed and wavelength as they approach the longitude of the 
jet velocity maximum, achieve their maximum amplitude a 
considerable distance downstream of the jet maximum, and 
eventually decay as the disturbances move out of the easterly 
jet regime. Furthermore, the high and low pressure centers 
of these transient wave disturbances occur on the wings of 
the central axis of the jet stream. Colton concluded that 
the transient disturbances are due to scale interactions 
involving short-term barotropic instability. 

In both Krishnamurti's observational and Colton's 
numerical studies, the zonal variation of the jet apparently 
has significant effects on the dynamic behavior of the 
transient disturbances. This study is an attempt to under- 
stand better such effects. This is a linear study, therefore, 
the results can only be suggestive of the complicated nonlinear 
behavior of the real atmosphere. However, since the ampli- 
tudes of synoptic waves are generally of smaller magnitude 


117 5 Opics than in middle and higher latitudes, a linear 








stability approach should be a reasonable approximation to 
the real atmosphere. 

This study is an extension of the linear barotropic 
instability theory developed over the years by various 
investigators. The extension of this study is the zonal 
variation of the basic flow. .Rayleigh (1890, 1913) developed 
the concept of hydrodynamic instability and Kuo (1949) 
extended this concept to a rotating atmosphere by including 
the beta term. Since then, many investigators have studied 
5277 bili oy of barotropic zonal flows. Kuo (1949, 1951), 
Lipps (1962), Yanai and Nitta (1968) investigated barotropic 
instability of parallel symmetric westerly jet flows. Nitta 
and Yanai (1969), Yamasaki and Wada (1972), and Kuo (1973) 
extended this theory to parallel symmetric easterly jet flows. 
Lipps (1965, 1970) and Kuo (1973) further extended the concept 
to parallel asymmetric barotropic zonal flows. In general, 
when the necessary criterion for linear barotropic instability 
is satisfied, normal mode solutions for barotropic unstable 
waves have certain properties: 1) they exist within a range of 
intermediate wavelengths, 2) they have a latitudinal tilt 
opposite to the shear of the barotropic zonal wind and this 
tilt creates momentum fluxes which redistribute the kinetic 
energy from the mean zonal flow to the disturbance, 3) their 
phase speeds in a westerly jet are less then the speed of the 


maximum westerly wind, and 4) their phase speed in an 





easterly jet may be greater than the maximum easterly wind 
(Pedlosky, 1964; Yamasaki and Wada, 1972). In this case, 
however, the growth rates are generally very small. 

Nitta and Yanai (1969) modified the concept that 
unstable solutions exist only within an intermediate range 
of wavelengths. For an easterly sinusoidal jet flow, they 
found a distinct short wave cut-off but no apparent long wave 
cut-off for instability. Yamasaki and Wada (1972) modified 
this by showing that the long wave cut-off is dependent on 
the strength and sharpness of this sinusoidal easterly jet. 
For very strong and sharp velocity profiles, the long wave 
cut-off approaches infinity. For relatively weak but still 
unstable profiles, a finite long wave cut-off exists for 
this jet. Kuo (1973) referred to this long wave cut-off 
region as the modified Rossby regime in his numerical study 
for an easterly Bickley jet. 

All these stability studies have used parallel baro- 
meopic flows, 1.€., there is no Zonal variation in the basic 
flow. Lorenz (1972), however, investigated primarily by 
analytical means the barotropic instability of a flow pattern 
which varies with longitude. The basic flow is a neutral 
Rossby wave superposed on a uniform westerly flow. Essen- 
tially, this flow depicts the progression of large scale 
waves embedded in a westerly current. Zonal flows of mid- 
latitudes are generally considered to be barotropically 


Stable, but Lorenz showed that a uniform zonal flow together 
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with superposed neutral Rossby waves may be unstable with 
1 CC но Tube: perturbations. Lorenz found that the 
growth rate of the perturbations is comparable to the growth 
rate of the errors of large numerical atmospheric models. 
Based on this, Lorenz suggested that barotropic instability 
may be partly responsible for the unpredictability of the 
real atmosphere. 

In addition to the regular or normal mode eigensolutions, 
there are "singular" or "continuum" mode solutions which 
have been discussed by Rayleigh (1913), Case (1960), Pedlosky 
(1964), and Yanai and Nitta (1968). These "singular" 
solutions correspond to continuous eigenvalues of phase 
velocity (c) which are equal to the basic flow, U(y), some- 
where in the zonal current. Since these "Singular" solutions 
are continuous, there is an infinite qübar of solutions. 
Case showed that these continuum modes are needed to form a 
complete set of solutions. He also showed that the disturb- 
ance formed by the sum of the continuum modes has a 
y-structure tilt in the same sense as the basic flow shear 
and that it usually decays as I/t or faster, where t repre- 
sents time. 

It is clear that in the real atmosphere, "mean flovs" 
vary both in space and time. They are neither purely 
baro roD ii nor purely baroclinic, and in including both 
5200 “010 a linear stability study is a very difficult task. 


Most of the investigators have examined the dynamic 





instability problem by studying each effect separately. 
This study, however, focuses only on barotropic instability. 
We hope that a better understanding of the kind of upper 
tropospheric waves studied by Krishnamurti (197la, 1971b) 
and Colton (1973) near the easterly jet can be achieved by 
studying the behavior of barotropic waves in a region of a 
variable mean easterly wind. We are especially interested 
in the situation where the waves move into and out of baro- 
tropically unstable regions. Therefore, the objective of 
this work is to determine the dynamic stability behavior, 
the structure and the energetics of barotropic waves propa- 
gating in a zonally varying mean wind. 

With a very simple analytical model, Pedlosky (1976) 
studied the finite amplitude dynamics in a zonally varying 
baroclinic current. The mean flow of his model changes 
abruptly from a weakly unstable regime with a constant growth 
rate to a stable regime downstream. He found that disturb- 
ances may propagate into stable regions substantially 
undiminished, retaining a considerable memory of its history 
in a locally unstable region. 

The mean flow in our barotropic study, however, has a 
full downstream variation such that the stability of the 
mean flow ranges from locally strong instability to locally 


dong stability. 
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Our approach is as follows: firstly, we develop a 
numerical model based on the linear non-divergent barotropic 
vorticity equation on a beta-plane incorporating the dynamics 
of transient barotropic waves in a region of variable mean 
wind. The mean zonal wind is an easterly hyperbolic secant- 
squarred (Bickley).jet and the mean meridional wind is derived 
in such a way that the mean flow is non-divergent. Secondly, 
using this numerical model, we conduct selected case studies 
to determine fundamental dynamic stability and energetic 
properties of these barotropic waves. Thirdly, we compare 
the results obtained from the numerical model with the 
parallel flow theory. This is done by constructing a simple 
mechanistic analytical model which incorporates the local 
stability concept of the parallel flow theory. This compari- 
son gives us further insights on the effect of the variable 
mean wind. 

The domain of the numerical model consists of an open 
channel with rigid walls at the north and south boundaries. 
The basic flow is a slowly varying easterly Bickley jet. 

A periodic forcing is applied on the inflow or eastern 
boundary such that periodic perturbations are generated from 
this boundary into the channel flow which represents waves 
moving into the region from the east. The western boundary 
condition allows the waves to move out of the channel. As 
the waves move through the region, they grow or decay in 
relation to the local stability properties of the mean flow, 


whereas at each point the fields vary periodically. 
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Il. NUMERICAL MODEL 


Ar GENERAL FORMULATION 
The governing equations are the barotropic vorticity 


equation using a beta (8) plane approximation: 


d 
d 


Qə 
Уҹ 
Уҹ 


du 


d 
+ .. + gv = Q - Dec 52.10) 


ct 
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and the non-divergent continuity equation 


ou öv _ 
am Ü : (2.2) 
where 
500”. Su 
ə öy (eS) 


Here ß is the north-south gradient of the earth's vorticity 
and is given by its value at 10° latitude. Q is a forcing 
function representing non-barotropic effects which is 


required to maintain the vorticity field, D, is a frictional 


m 
coefficient. Since the flow is two-dimensional and non- 


divergent, a streamfunction (y) is defined such that 


Q2 
= 


3v - 
u = - mu V (2.4) 


Qo 
X 


thus 


BD, (2.5) 
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Here u and v are the velocities of the flow in the x and y 
direction, respectively. Equation (2.1) now becomes 


2 2 : 
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The governing equations are linearized by separating 


quantities into the following form: 
2:7?) El” y ya E). | (8.3) 


Here the bars are the basic state variables and the primes 
are the perturbation quantities. The mean streamfunction 
v(x,y) satisfies the mean vorticity equation when a non- 
barotropic source term is added. The basic state equation 
is therefore 


9779” 90 377 
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d 
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3$ 28. - 
m DESC Q D ¿Y dë ла (2.8) 


Ke 


which is non-linear. Ü is needed to maintain the basic 
vorticity field in steady state. The resultant linearized 


perturbation equation is: 
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The advective terms in (2.9) can be written in Jacobian form 


such that 
99” - 2 zı 2 1 2 
ao * - Jiy,v V')- J(y sv v) - T SEAN 


where (2.10) is a Poisson equation with the tendency of the 
perturbation stream-function as the dependent variable. The 


right side of (2.10) is the forcing or non-homogeneous part 
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of the Poisson problem for ay'/3t. The north-south boundaries 


are rigid walls where 


ys Ü,” (2) 
and 
477 
ay * 0%, 12212) 


The inflow or eastern boundary conditions are specified such 


that 


$(0,y,t) = Aly) sin(ut) + 8(y) cos(ut), iS, 


and 
397 (0,y,t) - A(y)acos(at) - s(y)asin(at), (2.14) 


where w is the specified frequency of the forcing. A(y) 

and B(y) are the y-structure coefficients of the forcing 
which are described in Chapter IV. The purpose of this 
forcing is to introduce waves into the region from the east. 


The other condition on the inflow boundary is 


c'(O,y,t) = rept = - key e OH al) 


öy 
This expression follows because it is expected that y' will 
have a spatial variation of the form u The determination 
of k is discussed later. This boundary forcing is expected 
to lead to periodic wave disturbances throughout the domain 
after a certain time period of numerical integration. It 


turns out that the outflow boundary conditions are crucial 
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2::::11nT1fg this periodic state, because the wave disturb- 
ances must be able to propagate freely out of this boundary. 
Therefore, a Sommerfeld radiation condition 15 s 7) 
157 Dosneche tendency and the vorticity is used to approxi- 


mate this mechanism on the outflow boundary: 


Bea PS a (aye 

ot Ke So ox Bee ), and (2.16) 
a oz” 

uu ss "a 3x ^" 1501 


Here С̧и 15 а specified constant phase velocity. Pearson 
O 


(1974) showed that the Sommerfeld radiation condition is a 
consistent boundary condition for numerical models of initial 
value systems admitting dispersive waves. The fundamental 
problem with most outflow boundary conditions is the reflec- 
tion of incident waves from a boundary back into the interior 
region. This is usually disastrous for numerical models 


admitting dispersive waves. Pearson showed that, if d is 
O 


chosen judiciously, the longwaves move smoothly through the 
boundary. For short wavelengths, however, an area with a 
large coefficient of viscosity near the boundary is often 
required to help control the reflection problem. According 
to Pearson, the amount of damping is proportional to the 
wavenumber. A Matsuno finite difference scheme is used in 
our numerical model. This scheme has a tendency to damp the 
short waves thus we need only to use a relatively small 


Bere rent of friction to control wave reflection at the 
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gü í ow boundary. A friction coefficient of H = ERE” 
sec”! is used which is equivalent to an e-folding decay time 


of approximately 7.7 days. 


Bi. ZONALLY VARYING BASIC FLOW 
The basic velocity field is an easterly Bickley jet 


defined by 


ulx,y) = - U(x) sech“ (227) p m. 


Here d(x) is a characteristic length scale of the jet and is 


= 
| 
Q? 
* 


(20 


Qə 
< 


related to the half width d(x) (Kuo, 1973) by 
He) = 21.76 dix). (2.19) 


U(x) is the velocity of the jet at y=0 and U, Ís a constant 
velocity. The basic flow streamfunction (y) is specified 
Bomber constant at both the southern (y=-D) and northern (y=D) 


boundaries as follows: 


*(-0 =" 
Bb , , where (2.20) 
y(D) = C, | 
y 7 (2221) 
B(xsy) = U(x) | sech”(gyiylay s üç fay £ 9(-0), 
-D -D 
and 
D D 
DI = U(x) Ј sech^ (374) dy + U | dy 
Gap E 
EDO EC (2.22) 
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por a symmetric jet profile about y=0, C = (2. Equations 


(2.21) and (2.22) are conveniently simplified: 


Миу) Ud) (tanh (arxy) + tanh (gd) 
(2.23) 


+ Uy + UD + Ӱ(-0), 
and 
V(0) - 20(x)4(x) tanh (yö) s 20,0 8 9(-0) (2.24) 
From (2.23). the jet velocity at y=0 is given by: 
Y pi - 20,0 - 9(-0) 1 D 
55- ил ил J GO En “qu 5705) 


Therefore, if d(x) varies slowly in x, then, so does the 
basic flow. The x-variation for the characteristic length 


scale is given by 


ә 


7757 : 
850 km + 350 km 4cos[2* — | e X 2 XQ 
L = 


d TONES 
e 
О 





1200 km ə 


|^ 


Here 7 1s the longitude where the x-variation of the cosine 
function starts and L is the wavelength of this variation. 
Figures 1-3 show the basic fields of streamfunction 
fier, ZONnal velocity U(x,y) and vorticity c(x,y), 
respectively. Here L is set to 43,000 km. The domain is 
40,125 km long and 4,000 km wide (D = 2,000 km). At y=0, 


1 


the basic zonal velocity is -15.51 m sec " at the inflow 


boundary (x = 14,625 km) increasing slowly to a maximum value 
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of -30 m sec at x=0. From this longitude the central 


velocity slowly decreases downstream to a minimum value of 


I at x = -21,375 km. Between this longitude and 


-13.4 m sec 
the outflow boundary (x = -25,500 km), the basic flow is 

parallel. From Eq. (2.26), we note that the characteristic 
length scale, d(x), varies between 500 km at u(x,o0)nax and 


1200 km at u(x,o), i, 
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eO AS T AENT ELT OF BASIC FLOW 


The local growth rate of the variable mean flow of the 
numerical model is first determined by a parallel flow (i.e., 
no x-variation) numerical model of Williams et al. (1971) in 
order to gain some insight on the stability characteristics 
of the mean flow. This model is hereafter referenced as the 
parallel flow model. By setting au/3x = O in (2.9) of the 
numerical model, the governing equation of this model is 

арх ар Јана la 208] 29^ Lg. (3.1) 

at 3X Í aye 9 X 

The symbols of (3.1) have the same meaning as the numerical 


model of Chapter II, except that the basic flow is given by 


u(y) = - U sech (X) ; (3.2) 


where U is a specified constant that scales the magnitude of 
the central velocity of the Bickley jet (y=0). The character- 
istic length d is also a specified constant. Assuming that 
all perturbation quantities are periodic in x, Eq. (3.1) is 
finite Fourier transformed in x with wavenumber k, and is 
solved as an initial value problem. This approach gives the 
phase speed, growth rate, and the wave structure of the most 
unstable mode. It is convenient to write y' in the following 


form: 


u'lx,yst) = Alyst) cos(kx) + Biy,t) sin(kx). 6575 


A9 





Equation (3.3) is substituted into (3.1), and the coefficients 
of coskx and sinkx are separately set equal to zero. The 


205 term gives: 


2 E 
3 253A д 2 ð u 
Paes lage = rl er kt — B 
öy 7 öy öy 
5 | (3.4) 
“иә (Јл 
"zı 
öy 
and the sinkx term gives: 
2 2- 
| 2498 - 2 g ð u 
[37 - k$122 » KÉ (3 - k”) + (8 - Sa 
öy at au "E 
2 2 
L rer 50 
ду 
The boundary conditions are 
mu aS bl C. Dat) =O. (3.6) 


Equations (3.4) and (3.5) are written in finite difference 
form such that the second derivative with respect to y of a 


typical variable A is approximated as 


2 ün - 2A. + A. 
(251 z EE | (3.7) 
ду Ay 


J 
where j is the y-grid index and Ay is the distance between 
grid points. The model has 32 grid intervals and Ay is set 
equal to 125 km (the width of the channel is 4000 km). This 
y-grid structure is the same one used in the complete 
numerical model. Centered time differences are used for all 


quantities except those involving friction. The time step 
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(At) is set to 1 hour. The friction terms are evaluated at 
the previous time step in order to insure linear computa- 
tional stability (Haltiner, 1971). The integration begins 
with a forward time step. The boundary conditions (Eq. 3.6) 


become 


A = B = 0 at y = +D. (3.8) 
Equations (3.4) and (3.5) are solved for the tendencies by 
the exact method of Richtmyer (1967). These equations can be 


solved numerically as a function of 8, D u(y), and k for 


e: 
any initial conditions. These equations are integrated until 
their solution becomes exponential in behavior. The basic 
flow u(y) is specified for a number of selected longitudes 

of the variable mean flow of Fig. 2. The wavenumber k is 

also specified for each integration. 

In general, these equations have a set of discrete normal 
mode solutions as well as a continuous spectrum of solutions 
(Case, 1960, Pedlosky, 1964 and Yanai and Nitta, 1968). Only 
the normal mode solutions can give significant growth and the 
most unstable mode will dominate after a sufficient period 
of time. 

This parallel flow model is numerically integrated to 
150 days. The initial disturbance amplitude has a north- 
south structure with no tilt. The eigensolution obtained 
includes wave structure of the most unstable mode, and its 
2200000 ce and pase speed. By selectively testing different 


wavelengths, the most unstable wavelength is determined. 


3) 





Figures 4-7 illustrate the eigensolutions of the most 
unstable discrete mode for the selected jet profiles. Since 
the jet flow and the wave disturbance are both symmetric 
about y=0, only the lower half of the y domain is shown. 

It is interesting to note that the eigensolution amplitude 
Ras three maxima, one at y=0 and one on each wing (only lower 
wing is shown) of the jet approximately 600-800 km from y=0. 
It is apparent that when the jet is relatively sharp (i.e., 

d < 700 km), the central maximum predominates. When the jet 
is relatively smooth and broad (i.e., d = 800 km), the 
maximum on the wings predominates. Figure 8 shows the 

growth rate (n) corresponding to the most unstable wave- 
lengths (L) as a function of x based on the parallel flow 
model. In a local stability sense, the most unstable 
wavelengths range between 3650 km at x=0 where the jet has 
its maximum velocity and 4600 km near the inflow and outflow 
regions of the jet regime. We observe that locally the 
Gpeacest instability based on both tilt and growth rate is 
indicated where the jet achieves its maximum central velocity 
(u = -30 m sec"). This parallel flow model can only give 
the most unstable discrete mode; it can not depict any 
dynamic damping. This is reflected in Fig. 7(b) which 
corresponds to x = -22,500 km, outside the unstable region. 
Here the solution exhibits essentially no tilt and the growth 
rate asymptotically approaches the linear frictional damping 


25 unis is further discussed in Chapters VI and VII. 
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These results are consistent with those calculated by 
Kuo (1973) in a numerical study of a parallel Bickley jet 
and, in a local stability sense, may be viewed as a first 
approximation of the behavior of moving waves within a 


zonally-varying mean wind. 
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AFINE DIEFERENCE EQUATIONS 


The governing Eq. (2.9) of the complete numerical model 


is written with the advective terms in Jacobian form: 


397 20 
ei Bax = Dev Y 5 (4.1) 


an. ga ab ə 


for any two scalar quantities a and b. The Arakawa (1966) 
finite difrerence approximation for the Jacobian is used. 
This scheme conserves both the mean Square vorticity and the 
mean kinetic energy. 

We define the tendency of the disturbance streamfunction 
by 


L ду)" 
R 


> (4.3) 


where superscript t is the current time, and subscripts 1 
ana] refer to the x and y grid points, respectively. In the 
ANO ina finite difference equations, At is the time step 

and is set equal to one hour. Ax and Ay are the x and y 

grid point intervals, respectively. Equation (4.1) is written 


5007 16 İatsuno or Euler-backward finite difference scheme: 


2_t k t AO Ú ı 2- 
Y Tij ES eee Y ] ga. y] 
; (4.4) 
t 
y! “ui t-At 
es ihi. ” 
| 224 "adını 





177 Dee DEE (4.5) 


2 * * * 


me Ge EE 
V TaJ e -ell: 19.7 Y 1 all, Di m y] 


x * 
2. "E 
z Pina ins SE 
zı 224 |. с, ә 
15: = "с E 
vi; E es дә. (4.7) 


The finite difference form for the Laplacian of the tendency 


il given by: (4.8) 
t Acne t er t t 
grt. n aoc R5 qu E 
WE ac Ay” 


500 Oor the Jacobian terms, it is given by: 


= ü (y - y y (a = E ) + 
, j “AXaAy 177 LES i, j+] iE 
(4.9) 
ih T t Hr 
(Yi ja = Pa j-1) (5161.3 = 511,3) | 
(2 t 
z mül n int _ E 
ell, ; xay Си 7177 a eo) 
T, кс t 
i Yig, j EiT j+ mx ви S 
t ° 


cc MU 


: 4 t 
N Dear E = lp]. 
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@t 


: (t l | 
ell `" mm ERS ae E (4.11) 


1 t eg 
A E 


) + 


y T T E 

= тз 58,46 " Dat, 241" 
Ue T m 

аз 2 Мәлә) 


L qe | qe | = 


1:1 e? 19 i 
3 


(4.12) 


ell (9 sj is expanded in the same manner as for 
i,j 
Ë 


I [$,7^9'] 


1,J 
The Poisson Eq. (4.1) is solved for the tendency with 
a direct method developed by Sweet (1971). This direct method 
solves a finite difference approximation to Poisson's equation 
on a rectangular domain with Dirichlet boundary conditions. 
The finite difference form of the boundary conditions are as 


follows: 


a) Northern (j=J) and southern (j=0) boundaries: 


Vi - 0: i=0,I for j=0 and j=J, 25 
57:- i 
jo * 3 1=0,1, (4.14) 
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EE (4.15) 


b ) Ma on orea tern boundary 11, 1-0, 1) 


-> 


Up = A; Sin(weat) + B cos(u£At), (416) 
EN = A,ucos(wLat) - B. usinfuzAt), 1:17) 
and 
t t Ç 
uu 20205: 
meee tee 111 l, 1.3-1 (4.18) 
zs JJ Ay” 


where w is the specified frequency and k is a specified 


I ,J 
wavenumber structure. Time is discretized by: 


w ee 0 € 01.Z.3.4 ... (4.19) 


e and B. ANSeS ə rue ur coesticients of tne forcing 
(4.14). These coefficients were determined from the eigen- 
solution of the parallel flow model. The velocity profile 
for determining A. and E of the inflow boundary condition 
was selected from the inflow region of the numerical model. 


iki profile is described in Chapter V. 


с) Outflow or western boundary (i=0, j=0, J) 


E t 
o o. Di, EN (4.20) 
0,j I+c_ At/AX Ipaa NEUE ә : 
r r 
O O 
1 -At 
it ни E 2M AUS 
7752050 Ee HO o i 
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t ICAL 


* 0,3 ^ "0,j ди. (4.22) 


t 
The tendency To : 


2 


717721) and, Co, (4.21) were derived from 
the analytical Sommerfeld radiation condition (2.16) and 
(2.17), respectively. Since (2.16) is not in the usual form 
wene Dirichlet boundary condition, the direct Poisson 
Semetion solver, POISDD, had to be modified. Details of this 


modification are given in the Appendix. 





ДЕ АР СК ЕММИ 


Тһе investigation of the numerical model of the variable 
mean wind comprises two experiments. Experiment I is the 
principal one and will be discussed in detail in this chapter. 
In Experiment II a constant velocity is added to the mean 
1277 nd its results will be discussed in Chapter VII. 

The forecast equation 1s integrated in time from an 
BB state of ç ” “ 0. Tne periodic forcing on the eastern 
boundary causes the interior streamfunction to grow and the 
meor aton is continued until the time variation is periodic 
everywhere with the forcing frequency. By day 70, the model 
has achieved the fully periodic state and the wave packet 
envelope has become quasi-stationary. The forcing frequency 
is varied until the value which gives a maximum perturbation 
amea tude is found. 


The following values are used in Experiment I: 


2D = 4000 km, xp-x, = 40125 km, ü(0,0) = - 30 m sec^!, 


U 0, L = 43000 km, Ax = 375 km, Ay = 125 km. 


0 
For these parameters the maximum response occurs for a 
coe period of 3.25 дауз It requires several] complete 
integrations to refine the inflow and outflow boundary 
conditions for the specified forcing period. Wavenumber (k) 
for the inflow vorticity boundary condition (4.13) was 
determined by observing the predominant wavelengths near the 


inflow region. The radiation phase velocity ES ) in the. 
0 
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outflow boundary condition (4.15) was similarly determined 
by observing the predominant wavelengths near the outflow 
region. Since the domain is periodic, the phase velocity is 
approximated by the simple relation: 

| E kc ; (511 
Here o is the local frequency which,in the fully periodic 
istate,1s equal to the forcing frequency, and K. is the wave- 
number measured near the outflow region. Thus, for Experiment 


I, a phase velocity ER ) of -13 m 5 1s used for the 
0 


outflow Sommerfeld radiation boundary condition. 


A. RESULTS 

The solution for the Case I experiment becomes fully 
periodic after 70 days. The »'.field at t=70 days is shown 
in Fig. 9. An entire train of barotropic waves actually 
exists throughout the lengtn of the channel, but the waves 
upstream of x = -4,500 km and in the outflow region are not 
shown in Fig. 9 because of their relatively small amplitude. 
ie maximum Wave amplitude occurs at x = -12,750 km and is 
2 and 4 orders of magnitude larger than at the jet maximum 
(x=0) and the inflow boundary (x = 14,625 km), respectively. 
Figure 10 shows the lower half of phase angle tilt 9*(x,y) 
of the wave disturbance y-structure for various values of 
longitude x. We note that the waves upstream of approximately 
x = -20,000 km tilt opposite to the mean wind shear, which is 


15757 7 barotropic instability. In fact the tilt near 
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the inflow boundary shows relatively weak instability and 

this tilt slowly increases, reaching a maximum at x=0 where 
the jet velocity is maximum. Further downstream, the tilt 
slowly decreases and eventually reverses near x = -20,000 km. 
This behavior is consistent with the growth rates for the 

most unstable wavelengths shown in Fig. 8 which were computed 
using the parallel flow model. Near the outflow boundary, 
however, the tilt is reversed. This indicates dynamic 
2150771167 or a Tlov of energy from the disturbance back to 

the mean flow. Therefore, we observe that the tilt of the 
wave disturbance qualitatively adjusts to the local stability 
of the mean flow (see Figs. 4-7). We recall that the parallel 
flow model can only solve for the most unstable discrete mode, 
220) dynamic dambing i$ not indicated in Fig. 7(b) for the 
guet low region. 

Figure 11 shows the envelope of the wave packet, 
evaluated at y = -/50 km, where the disturbance amplitude is 
large. This envelope <w'(x)> is obtained by recording the 
maximum and minimum y' values that occur at each longitude 
over a period of ten days after the solution becomes fully 
periodic. If a larger time interval were chosen, the envelope 
would not change. Note that the maximum amplitude occurs in 
the area where the local growth rate becomes zero (see Fig. 8). 
The smoothness of the fields in Figs. 9 and 11 indicate that 


the radiation outflow boundary condition is working properly. 
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Тһе ф" fields for days 83 through 86 are shown in Figs. 
12 through 15, respectively, to illustrate a time sequence 
Eu new banrotropic wave train. A periodicity of 3.25 days 
can be determined by examining the time series at each grid 
point (not shown). We also note from this sequence that there 
are three maxima of y' in the wave structure. One maximum is 
at y=0 and another is on each wing of the jet near y = +750 km. 
The maximum on the wing clearly predominates. Figure 16 shows 
the variation of wavelength, L(x), for y=0 and y = +750 km 
corresponding to the latitudes of the observed three maxima 
of y'. The disturbance wavelength is about 4100 km initially 
near the inflow boundary. The wavelength near y = +/50 km 
is larger upstream and smaller downstream of x=0 than at y=0. 
At y = +750 km, the maximum wavelength of 5060 km occurs 
approximately 950 km upstream of x=0 while the minimum wave- 
length of 3600 km occurs nean the outflow boundary. At y=0, 
the maximum wavelength of 4900 km occurs about 950 km down- 
stream of x=0 while the minimum wavelength of about 4000 km 
occurs near x = -13,500 km with the wavelength increasing 
slowly further downstream to the outflow boundary. The 
range of the wavelength is between approximately 3600 km and 
5000 km, which is nearly the same range of values obtained 
for the most unstable wavelengths using the parallel flow 
model in Chapter III. Figure 17 shows the disturbance phase 
speeds, 5P: for latitudes y=0 and y = +750 km where 


5 is obtained from 
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Си = L(x)w (5.2) 


Since the entire domain is periodic with the forcing 
frequency, gx) has the same basic behavior as L(x) 
[Fig. 16] for the corresponding latitudes. Namely, upstream 
of x=0 the disturbance phase velocity is faster on the wings 
of the jet than at y=0, and downstream of x=0 it is slower. 
This behavior of both L(x) and ae illustratesthe-til 
behavior of the waves as they move into and out of the baro- 
tropically unstable region of the variable mean wind. More 
results are presented in the next chapter. 

The parallel flow model is also employed to determine 
a reasonably good forcing function that would generate 
pertodic wave disturbances from the eastern boundary. This 
is accomplished as follows: 

Equation (2.14) defines the analytical form of the 
periodic forcing which is applied as an inflow boundary 


condition in the complete numerical model, i.e., 


a (0,y,t) = Aly)ucos(wt) - B(y)usin(ot) 


Equation (4.17) describes the finite difference approximation 

to this forcing. The coefficients A(y) and B(y) of (2.14) 

are obtained from the eigensolution determined using the 

parallel flow model for the following Bickley jet profile: 
(y-0) - .16 m sec^! 


u 
d = 962 km. : 
L = 4,600 km 


5 





This profile is representative of the mean flow near the 
inflow region. Figure 18 shows the y-structure of these 
coefficients. 

Two other forcing functions, which are independently 


constructed from the parallel flow model, are also tested. 


These functionsare: 


=- O cos (52) sin(wt), 55 
and 

EK T 8 

SE (0,y,t) = - o (cos(—>Z)) sin(ut), (5.4) 
where D is the half width of the channel (D = 2000 km). 


The behavior of the wave disturbances for all three 
forcing functions are essentially the same except in the 


vicinity of the inflow boundary. The parallel flow solutions 


for the most unstable wavelength show a weak decay rate in 


the inflow region (see Fig. 8). This is because the linear 


friction decay rate is 
growth rate. We found 


the inflow boundary is 


slightly larger than the weak local 
that this weak decay behavior near 


achieved in the numerical model using 


socia function (2.14). Use of the other two forcing 


functions, (5.3) and (5.4), realizes a weak secondary <y > „x 


near the inflow boundary. Therefore, we found through use 
of the parallel flow model that the first forcing function 
provided the better inflow boundary condition than the other 
two, although no substantial difference occurs even if the 


other forcing functions are used. 
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Fig. 18. The coefficients A(y) and B(y) 
of the eastern boundary forcing 
function used in the complete 
numerical model. 
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B. ENERGETICS 

To understand the dynamics of barotropic instability, 
it is helpful to derive the disturbance kinetic energy 
equation. A finite difference approximation form of this 
equation is used to compute the energy budget of the waves. 


The linearized momentum equations are: 


Sup" aud au E I 20 


m ә фај "ју T 
| 55 
- syv' = - — 22 - oyu’, 
р 
and 
3v' SERA 3V 2 „av! ER 
BB su m aem M s 
50) 
P yu ie m Se - Dev’ 
р 


25 /71ltipiyinu (5.5) by u” and (5,6) by v' and adding the 
resultant equations, we obtain the disturbance kinetic energy 


equation: 


ok" =- ð -3k ' =. s 7:07 ! non e 1 (9v + 
ә — e DE e u' v Jy E 
, | 7577. 
ити Ји mın bi 
əy - aX - öy f 
p p 
where 
T Ta 
59 x. — 


The disturbance fields are considered to have a time 
dependence of the form sinlat + e(x,y)]. Thus, when 156) 
is averaged over one period and over one latitudinal domain 
(-D<y<D), the following energy balance equation is obtained: 
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1 2 3 4 5 
EN “аб "ә 7 “əy " 
3X 3x Der 


5: 
P 
Term ] is the advection of k' by the mean flow (u). Terms 

2, 3, 4 and 5 are the barotropic exchange terms with the 

mean flow. Term 6 represents the work done by the pressure 
mud Term 7 is the dissipation term and represents an 
energy sink. With the exception of term 6, all the terms in 
5777 are determinable from the J and w' fields. Thus, the 
pressure work term is determined as a residual in Eq. (5.9). 
The energy equation is evaluated after the complete numerical 
model achieves the fully periodic state. 

Figure 19 shows the energy balance as function of 
longitude. Clearly, the maximum value of each of the terms 
in the energy equation occurs downstream of x=0. Curve (a) 
represents the Reynold's stress term, <(-u'v'au/3y)> , where 
the maximum value occurs approximately at xX = -10,975 km. 
Curve (b) represents the advection term, <(-u3k'/3x)>, of 
k' by the zonal current. Curve (c) represents the linear 
dissipation term <(2D.k')>. We observe that the largest 
dissipation takes place near the maximum amplitude of the 


wave packet envelope, <w'(x)> Curve (d) represents the 


max” 


pressure vork алт ¿(A Q It also has a maximum 
, p 


near <w'(x)> The remaining barotropic terms are not 


max” 
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included in Fig. 18 because they are 2 to 3 orders of 
magnitude smaller than «(- u'v'3u/3y)». 


Recall that for Experiment I, <y'(x) occurs approxi- 


”max 
mately at x = -12,750 km. We observe that the advection of 
k' by ü plays a dominant role in the energy balance. Once 
the fully periodic state is achieved, the generation of 
disturbance energy by barotropic instability essentially 


occurs in the region between the <y'> (x = -12750 km) and 


max 
the jet maximum (x=0). It is in this region that the 
strongest positive correlations occur between u'v' and 
-3u/3y in the Reynold's term. This is consistent with the 
observed behavior of the waves (see Figs. 10 and 11) since 
this strong correlation of u'v' is dependent on maximum 
amplitude of y'squarred and the maximum tilt of the wave 
structure with the shear of the mean flow -au/ay. It is 
also from this source region that u transports a large part 
Mics generated k' downstream to a point centered near 
2 -15,.000 km. 

Strongest dissipation of k' occurs in the region centered 


near «y'(x) It is in this region that the strongest 


"max: 
spatial gradients of oi occur. The largest values of pressure 
work also occur in this region where the pressure work term 

is computed as a residual. We note that its location of 


occurrence is consistent with the large amplitudes and spatial 


gradients of oi of the region. We also note that о ла. 
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vanishes at <y'(x)> oy Since k' is maximum at this longitude. 
The disturbance energy balance which is shown in Fig. 19 
graphically illustrates the basic energetics involved in 
supporting the large amplitudes and spatial gradients of y' 

a considerable distance downstream of the jet maximum. 


EMI uv. che-Reynold's stress term, «pev a reene 


dominant source for disturbance kinetic energy. 
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WEY SIMPLE MECHANISTIC ANALYTICAL MODEL 


A. FORMULATION 

In this chapter we develop a simple mechanistic analyti- 
cal model which uses the locally determined parallel flow 
solutions. The results of this simple model can be compared 
with those of the complete numerical model to obtain some 
measure of the adjustment of the transient barotropic waves 
to the local stability of the variable mean wind. 

The following equation allows for propagation and growth 


or decay: 


(x) 93. 2 n(x)y' , (6.1) 


ow! 
at AC ox 


r 


where y' represents the disturbance streamfunction. c. (x) 
the local phase velocity and n(x) is the local .growth rate. 
We observe that if С̧и and n are independent of x, then 
21 7136 exact and, in generali, should give a reasonable 
approximation to the downstream variation of y'. 
Let us consider the case where y' has a solution of the 

form: 

y! s F'(x)e!"*, (6.2) 
where w is a specified frequency. After substituting (6.2) 
into (6.1), we obtain 

iuF' + c (x) SE = и(х)Е". (6.3) 


r A 
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Rearranging (6.3), ve get 


Bl dx E : (6.4) 


X 
à n(x) - i 
gn ет = f > doc Lone 775 
Xo 
" | | y (5.6) 
1 E 1 = m n(x) 
F (x) C | 1 J Sat 5” E f c (x) dxe 
X X BR 
0 0 
Some fundamental insight can be gleaned from (6.6). Here 


the amplitude of w' must be specified at some initial point, 
X» which could be the inflow point. The local wavenumber 
is о/с̧, and the spatial growth rate is n/c.: These are the 
real and imaginary parts of the wavenumber, if we were to 
write p'=A exp (ikx-iwt). 

In this simple analytical model, we are interested in 
the growth and decay of the wave packet. If the oscillatory 
behavior of Eq. (6.6) is dropped, the following equation 


expresses the envelope of the wave packet: 
m... 
m nx 
“50. en f —” (6.7 
"o 


Equation (6.7) gives the exponential growth/decay behavior 
of a wave disturbance as it travels within its wave packet 
envelope. This behavior is determined by the integral 


effects of the local stability properties of the mean flow. 
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D. APPLICATIONS 

Figure 20 compares the envelope for w'(x) from the 
complete numerical solution (Fig. 11) with the F(x) given 
ЫБ, (6.7) Тог two values of xg. In this figure the 
amplitude is plotted against upstream distance from the 
initial point E The analytical model uses the locally 
determined parallel flow growth rates (n) and phase speeds 
Le, based on the wavelengths measured from the complete 
momerical solutions. Therefore, Eq. (6.7) is applied at 
1507 5146 y = + 750 km where the wavelengths are measured. 
The curve for F(x) can be shifted up or down by changing 
x? but its shape does not change. In Figure 20 all the 
curves have a maximum amplitude at x= - 12,750 km where 
the local qrowth rate is zero; however, the simple model] 
has a lower maximum growth rate and slower damping rates 
than the complete numerical model. 

The lower portion of Fig. 2] contains the wavelength 
measured from the numerical solution at y = - 750 km. In 
126 Oper portion of the diagram are the phase velocities, 
C. from the parallel flow model and a from the complete 
numerical model for the same latitude. Here C, is computed 
with the use of the wavelength: shown. The two phase speeds 
have similar behavior although CT is shifted slightly 
upstream and has larger variations. We also note that oe 
is generally smaller than c, except for the region 


S500 km < x - 6750 km. 
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Similarly, Fig. 22 contains the wavelengths and phase 
speeds for latitude y=0. Again, the two phase speeds have 
similar behavior, but at this latitude ce is shifted 
slightly downstream and the differences between phase speeds 
is generally less than at y = -750 km. Again we note that 
— is generally smaller than c, except for the region: 
-2250 km « x « 2250 km. This indicates that both the wave- 
length L(x) and the phase velocity с̧*(х) between latitudes 
y=0 and y = -750 km are maximum where the mean flow is 
strongest, 

Figure 23 contains the spatial growth rates, m and m*, 
from the parallel flow model and from the complete numerical 
model, respectively. For the parallel flow model, we obtain 


m= = n/c, : (6.8) 


The value of m* is computed directly from the envelope. We 
note in Fig. 23 that m* has a larger maximum than m and the 
maximum is shifted slightly downstream from the jet maximum. 
Both curves pass through 0 at x = -12,750 km which is the 
maximum for both wave packet envelopes. However, m* 

shows a much larger damping in the outflow region which can 
also be seen in Fig. 20. In fact the parallel flow solution 
damps at the rate given by the frictional coefficient 
divided by the phase speed. This is expected in the baro- 
tropically stable outflow region when using the parallel 
flow model because this model can only give the eigen- 


solution of the most unstable discrete mode. 
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The solutions observed in the numerical model (see 
Fig. 10) toward the outflow boundary tilt in the same sense 
as the mean wind shear which gives dynamic damping. This 
tilt behavior is also indicated by the divergence of phase 
speed in the outflow region between y=0 and y = -750 km 
(see Figs. 21 and 22). This dynamic damping appears to 
be due to continuous spectrum solutions (Case, 1960; Yanai 
and Nitta, 1968) which are not included in the normal mode 
solutions that are employed in the simple integral 
[Eq. (6.6)]. In fact we observe that, in the outflow region, 
the dynamic damping is 1/t or faster, increasing near the 
outflow boundary. In Fig. 23 the m* curve is skewed slightly 
to the left with respect to m and the jet maximum. This 
could be the result of the tilt structure in the wave lagging 
spatially behind the expected value from the local stability 
conditions. This would give a smaller growth rate on the 
upwind side of the jet maximum and a larger growth rate on 
the downwind side. This effect is indicated by the asymmetry 
of the m* curve with respect to the m curve in the unstable 
region of the mean wind. However, this is not the only effect 
involved because the most striking feature of Fig. 23 is 
the fact that the maximum value of m* is significantly 
larger than the maximum value of m. We have already shown 
enat the most important disturbance energy production term 
(see Fig. 18) in the energy computations with the numerical 


Solutions is proportional to 


5 





7707? 01) 


This term is the only source term for the parallel flow 
model, and it depends on the phase tilt in the disturbance 
field. Figure 24 compares the phase tilt of the numerical 
model solution, 86”, and that of the most unstable vavelength 
of the parallel flow model, 0009 at x = 0, x = 3750 km, and 
x = 13125 km. It is evident that at the jet maximum the 
5313 from the numerical model is Significantly larger than 
no ne tilt from the parallel flow model. This is also 
evident but to a lesser degree at x = 3750 km. These results 
are consistent with the larger growth rates for the variable 
HB 100) hws stronger tilt is also indicated by the 
comparison of the Ee and c, curves in the region near and 
slightly upstream of the jet maximum for latitudes y=0 and 
y = -750 km (see Figs. 21 and 22). Ve note that at 
у = -750 km, C * is relatively larger than c. near the jet 
maximum. This difference is significantly reduced and 
shifted slightly downstream at y=0. He also note that near 
the inflow region (x = 13125 km), the phase tilt of the 
numerical solution lags slightly behind the parallel flow 
solution (see Fig. 24c). Apparently the downstream varia- 
tion of u augments the phase tilt which gives a larger 
moeow cn ate. 

In comparing the parallel flow theory and the complete 


numerical model, a source of error could be introduced 
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throuqh the zonal resolution. This is because in the former 
the zonal variation is specified by a spectral representa- 
sona while in the latter it is resolved by finite differences. 
Since the wavelength appearing in the complete numerical model 
varies between 10 and 13.5 Ax, such resolution normally gives 
a good approximation to the exact solution (Haltiner, 1971). 
117 this error should not be very significant. In order to 
confirm this we carried out an experiment in which a parallel 
mean flow is specified by the full two-dimensional model 
fintte differences of the complete numerical model. Тһе 
growth rate and phase speed calculated from the experiment 
agree very well with those obtained by the semi-spectral 
parallel flow model, the difference being only 1-21. Thus 

we may conclude that the differences between the numerical 
solution and the local solution of the parallel flow theory 
are genuine and are not due to the differences in the model 


resolution. 


Er PEEECTS OF LINEAR FRECTION 

Figure 23 shows that the local growth rate, determined 
by the parallel flow model, changes sign at X = -12,750 km 
and that it approaches the linear friction decay rate near 
the outflow region. We also note that the maximum aplitude 
of both the analytical and numerical model wave packet 
envelopes occurs at the lontitude where n(x) vanishes. 


Based on these results, <y'> is expected not only to 


max 
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decrease Tn magnitude but to occur further upstream if the 
linear frictional rate is increased. Likewise, the opposite 
effect would be expected if this friction were decreased. 


This hypothesis is tested by rerunning Experiment I using 


larger friction of De = 700 sec”! which is equiva- 


lent to e-folding decay time of 4.63 days. We found that 


<y > ay occurred approximately 1700 km further upstream and 


with an amplitude of an order of magnitude smaller (.4 vice 


5.9) than for the original experiment. 
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In Experiment II a constant advective speed of 
Дер -5 m 5 İs added to the mean flov, othervise the 
parameters are the same as those used in Experiment I. ‘We 
hope this modification will shed some light on the adjustment 
process of a faster moving wave chain. As in Experiment I, 
the analytical model also uses locally determined parallel 
flow growth rate (n) and the phase speed (c) based on the 
wavelength measured from the complete numerical solution. 
The local phase speed c is now related to the Doppler- 
shifted phase speed c by C = C + U The maximum w' response 
mnis experiment is obtained with a forcing period of 2.5 
days. Since the Doppler-shifted frequency may be written 
as R = k(c* - I = u -= SE this shorter period actually 
gives approximately the same Doppler-shifted frequency as 
that of Experiment I (where the period is 3.25 days with 
U = 0). Figure 25 shows the y' field for day 70. 

Based on the results of the simple analytical model 
[Eq. (6.8)], we expect that the constant advective velocity 
U. will significantly reduce the local spatial growth rate 


n 
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Figure 26 shows the spatial growth rates from the parallel 


flow model (m) and from the complete numerical model (m*). 
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Comparing this figure with Fig. 23, the reduction in both m 
and m* from those of Experiment I is clear. This reduction 
İs evidently due to the advective velocity which causes the 
disturbances of Experiment II to move approximately 

-5 m sec" faster than in Experiment I. Thus, the waves 

have less time to grow or damp in response to the local 
stability per unit distance traveled. The smaller spatial 
growth rate obviously affects the amplitude and spatial 
gradients of y'. Figure 27 shows the wave packet envelope 
<v'(x)> of the complete numerical solution for Experiment II. 
Figure 28 shows this envelope along with the analytical model 
envelope F(x) using either Xo = -375 km or ə 3750 km as 
Bio initial point (no discernible difference is noted between 
the use of these initial points). Both envelopes are, as 
expected, significantly smaller than those of Experiment I. 
In addition, there is less variation in the wavelength 

(3750 km < x < 4875 km) compared to Experiment I. This is 
also a result of the faster phase speed which causes the waves 
to have less time to adjust to the local u(x,y). 

Although the faster phase speed of the waves in 
Experiment II reduces the magnitude of growth and decay rates 
and therefore brings the spatial growth rate curves of the 
numerical model (m*) closer to that of the parallel flow 
model (m), we notice a stronger asymmetry in the m* curve 
with respect to x=0 in Fig. 26 compared to Experiment I 


(see Fig. 23). In addition to the maximum m* occurring 
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slightly further downstream of the maximum of m(x=0) than in 
Experiment I, the cross-over points of the m* and m curves 
exhibit more asymmetry in Experiment II than I. The first 
cross-over point occurs approximately at longitude x = 8500 km 
and x = 4000 km for Experiments I and II, respectively, while 
the second cross-over point occurs at approximately 
x = -12000 km where the spatial growth rates vanish. As in 
Experiment I, we interpret this asymmetry as being due to the 
tilt structure of the waves lagging behind the expected value 
mu local Stability condition. Ihe fact that this lag 
effect is enhanced by the faster phase speed of the waves in 
Experiment II can also be seen in Fig. 29 where the north- 
south phase tilt o*(x,y) of the wave disturbances for both 
experiments are compared at several longitudes. The unstable 
tilt in Experiment II is weaker upstream and stronger 
downstream of approximately x = -2000 km reflecting the 
stronger lag effect of the faster moving waves. x : 

Figure 30 compares the phase tilt of the numerical model 
solution, 8*, and that of the most unstable wavelength of 
the parallel flow model, 6, for two longitudes (x=0 and 
x = 3750 km). Comparing this figure with Fig. 24 of 
Experiment I, it is obvious that the apparent augmentation 
of the wave tilt due to the zonal variation of u(x,y) is now 
significantly reduced. This behavior may at least be 
partially explained by the lag effect caused by the advective 
current Un» which generally reduces the growth rate in the 


vicinity of the jet maximum. 
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The lower portion of Fig. 31 shows the wavelength 
measured from the numerical solution at y = -750 km. In 
the upper portion of the diagram are the phase velocities at 
y = -750 km, c* from the complete numerical model and e from 
the parallel flow model based on the corresponding numerical 
solution wavelength. The two phase speeds agree very closely 
upstream of the jet maximum while c* is smaller than c down- 
Stream of this longitude. Similagly, Fig. 32 contains the 
Wavelength and phase speeds for latitude y=0. At this 
latitude, the two phase speeds have similar behavior through- 
out the length of domain with c* being approximately 
Im sec! slower than c. Figure 33 shows the wavelength 
curves [L(x)] for latitudes y=0 and y = -750 km. At both 
latitudes the maximum wavelength occurs near the jet maximum, 
but the wavelength (and therefore phase velocity) at 
y = -750 km is larger upstream and smaller downstream of 
approximately x = -5000 km. As in Experiment I, these 
differences in c*(x) and L(x) between the two latitudes are 
consistent with the observed tilt behavior of the waves. 

Other than the foregoing mentioned differences, the 
basic behavior of the waves in Experiment II is similar to 
that of Experiment I. This includes the tilt of the waves 
247 - reversal Errom that Of growth to that#of damping 


downstream of SUE: a , the occurrence of SUME i a consider- 


ax 
able distance downstream of the jet maximum, the general 
shape of the wave packet envelope, and the longitudinal 


variation of the energy balance (Fig. 34). 
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VIII, SUMMARY AND CONCLUSIONS 


In this work we studied the behavior of waves moving 
into and out of a barotropically unstable mean wind field 
which varies both in x and y. This mean wind field roughly 
simulates the 200 mb easterly jet south of the Tibetan high 
during the Northern Hemisphere summer. A rectangular domain 
is used with a time periodic forcing on the inflow (east) 
boundary and a Sommerfeld radiation condition on the outflow 
(west) boundary. This allows the simulation of the propaga- 
tion of small amplitude waves through the easterly jet region. 
The vorticity equation on a beta plane is solved with the use 
of finite differences, and, when the boundary conditions are 
properly adjusted, the barotropic waves move smoothly across 
the region and out the western boundary. After a certain time 
period of numerical integration, the solution becomes 
periodic everywhere with the forcing frequency which is 
specified on the eastern boundary. As the waves move 
through the jet regime, they grow or decay spatially in 
response to the local stability of the mean flow, even 
though at each point the variation is purely harmonic. The 
numerical results are compared with a simple mechanistic 
analytical model which is developed using the local growth 
rate and phase velocity computed from the parallel flow 


theory. 
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Two experiments are carried out in this investigation. 
The parameters used in these experiments are the same, 
accept that in Experiment II a constant advective velocity 


. - ] 
of U, 5 m sec 


İs added to the mean flow in order to 
gain a better understanding of the adjustment process of a 
faster moving wave train. 

The basic behavior of the waves in both experiments is 
similar. This includes the following: 

1) The maximum amplitude of the waves occurs a 
considerable distance downstream from the most unstable part 
of the jet near where the local growth rate of the parallel 
flow theory vanishes. 

2) The tilt structure of the waves and its reversal 
ипи ат ок growth to that of damping occurs downstream of 
the maximum of the wave packet envelope in response to the 
local stability of the mean flow. 

ER The wavelength (and therefore phase velocity) 
increases slovly from the inflow boundary to a maximum value 
near the jet maximum and then decreases slowly downstream 
reaching a minimum value near the outflow region. 

4) Тһе Reynold's stress term <(u'v'au/ay)> is the 
dominant source for disturbance kinetic energy and the 
2116:61/€ term) -“(-üək /9x)? plays a Sigqnificant role 1n 
transporting a considerable part of this energy downstream 
from the source region. 

2) Large spatial damping occurs in the outflow region 


due to the presence of the continuous spectrum modes. 
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In Experiment I where the mean wind vanishes toward 
large y, the maximum spatia] growth rate in the numerical 
mode] is significantly larger than the spatial growth rate 
of the parallel flov model. Apparently the downstream 
variation of the mean wind augments the phase tilt of the 
waves thus causing an increased growth rate as compared with 
waves in a parallel flow. In Experiment II where the waves 
move westward faster due to the presence of a constant 


- 5m sec” 


advective zonal wind, the spatial growth rate 
computed by the numerical model is brought closer to that 
given by the parallel flow theory, but the asymmetry of the 
numerical growth rate curve with respect to the jet maximum 
is enhanced. The reduction of spatial growth and decay rates 
is mostly due to the faster phase speed which causes the 
waves to have less time to grow or damp per unit distance 
traveled. The enhanced asymmetry is apparently due to a lag 
effect wherein the adjustment of the waves is delayed 
waal ly. The Tatter effect is also reflected in the north- 
south phase structure of the waves and in the smaller varia- 
tion of the wavelength as the waves propagate downstream. 
The augmentation of the phase tilt by the downstream 
variation is also significantly reduced which may at least 
be partially explained by the lag effect. 

Recently Pedlosky (1976) used a simple analytical model 
to study the effect of the downstream variation of a baro- 


clinic current. The mean flow in his model changes abruptly 
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from a weakly unstable regime with a constant growth rate 

to a stable regime downstream. His model includes nonlinear 
effects but it does not produce strong dynamic damping due 

to the continuous spectrum modes. This is perhaps because 

he did not include the lateral shear of the mean flow. Thus, 
Pedlosky found that his disturbances remain substantially 
undiminished when they propagate from a locally unstable 
region to a locally stable region. 

The behavior of the wave disturbances in our study 
resemble certain aspects of the waves observed by 
Kraıshnmamurti (1971a, b) along the South Asia easterly jet 
during the Northern Hemisphere summer, and the waves that 
were simulated by Colton (1973) in a nonlinear model. The 
wave disturbances entering upstream of the easter]y jet 
regime grow and increase their speed and wavelength as they 
approach the jet maximum. They then achieve their maximum 
amplitude a considerable distance downstream of the jet 
maximum and subsequent]y decay as they leave the strong jet 
regime. 

The simple mechanistic analytical model developed in 
this study can crudely approximate the adjustment of the 
0 mg disturbance to the local stability of a variable mean 
flow within the dynamically unstable region. It is possible 
that this simple approach could be used to study the linear 


instability of other variable mean flows to wave disturbances. 
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Although this study employs linear equations, it is 
possible to estimate qualitatively the nonlinear effect of 
the waves on the mean flow. Since the Reynold's stress term 
is the dominant term in the energetics, most of the mean 
kinetic energy is removed by the waves on the downwind side 
of the jet maximum between where the wave amplitude is 
maximum and where the phase tilt of the waves is Strongest. 
Clearly these synoptic scale waves will affect the amplitude 
and phase of the quasi-stationary planetary waves which 
combine to form the easterly jet over South Asia. These 
effects can be studied by including nonlinear effects in 
the present model. It may also be interesting to extend the 
present model to a 3-dimensional model to study the effect 


pievertical shear and baroclinicity. 
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APPENDIX 
Subroutine POISDD 


Subroutine POISDD is a direct method solver of the 
finite difference approximation to Poisson's equation 


(Sweet, 1971): 


v^ u(x,y) 27:50 (1) 
for a<x<b 
C V qg 


and with Dirichlet boundary conditions 


u(x,c) = g(x) 

u mq) = g. (x) a<x<b (2) 
u(a,y) = g3ly) c<y<b 

ulb,y) = g,(y) 


We define a grid on the rectangle {(x,y): a<x<b, 


ee by selecting two Ppositivesinvegers M and N, such that 
b-a d-c 


AX = "NV and AY = mum (3) 
and where N must be a power of two, and the points 
x. = a + (i-1)Ax lin ə M 
(4) 
- = c+ (j-1)AY j  -. NI, 


where i and j refer to the x and y grid points, respectively. 
The Poisson equation that is solved by POISDD is given 


КР ТОЈ): 


E ı “550. 57... 2 d 
IA O = (OY 


57 





vith boundary conditions Eqs. (2.11), (2.14) ana (2.16). The 
finite difference form of these equations are given by Eqs. 
(4.10), (4.14) and (4.17). Here n is the finite differ- 
5055 representatTon for The tendency, 3y /at, and is the 
dependent variable in our Poisson equation. For our 
Numerical Model, M=106 and N=32; t is the current time, 
Bere tt = gAt 2 = 1,2.... and At in the time step. 

Our system of equations meets the requirements of 


subroutine POISDD, except for the outflow boundary condi- 


ones. (2.16). 


|o 


s RE). (6) 


Q2 


r 
O 


Q2 


X 


è 


This equation is the Sommerfeld radiation condition (Pearson, 
Ea Clearly, this is not the usual form of Dirichlet 
boundary condition. Therefore, in order to use POISDD in 

our numerical model, this subroutine had to be modified 

il orporatınq the radiation oútf low condition, This was 


10206 öy vriting Eq. (6) in Tinite difference form: 


t t-At t IU 
RS. 9 (7) 
At r AX ı 


where the finite difference symbols have their usual meaning 
in numerical weather prediction. Subscript i=1 refers to 
the outflow boundary, and subscript i=2 refers to the 

column of grid points next to this boundary. We observe 


that this scheme is backward in time with upwind space 
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differencing. 

that 

t-At 

NC AME 

eel 
0 


ХЫ) + “PS: 
e er Ax ] + 112: At) 


Equation (7) is conveniently rewritten such 


T t 


(8) 


O 


Since all the interior points of the domain must satisfy the 


non-homogeneous Poisson equation, we, then, use this constraint 


to couple the radiation boundary condition [Eq. 


(8)] with the 


interior points, as follows: 


L L 


T 215.1 + T 


AX 


The term inj 





t 


2A EN E 





(9) 
ZE à 
= elo j + T 


t 
PEN oe 

2.) 
Ay 


is the finite difference approximation to the 








right side or non-homogeneous part of Eq. (5) or Eq. (2.10). 
He multiply Eq (9) Dy E and we obtain 
t Í t 
I. 217 mə 
SEH 251 3 E L Ы t 
—. x... 
Ä (10) 
8 2 
= Тә jey 
where s = 2 
We now combine Eq. (8) and Eq. (10), and we get 
T,° + 2 + 2 - 5 SE r 
© 02,jtl S a 2: S 
Em At 
Es O (11) 
ә Ди 2 . 7 
“12.31 2 s[T+c, At/8x] 


О 


5 





By incorporating Eq. (11) in subroutine POISDD, we 
have, thus, coupled the Sommerfeld radiation boundary condi- 
tion to the interior domain. Also, by Eq. (8), we now have ` 
a complete set of Dirichlet boundary conditions. The finite 
Ae erence Egs. (4.10), (4.14) and (4.17), thus, form a 
linear system of equations of dimension (M-1) x (N-1) for 


the unknowns ə 2217777 17 7715 7? > Cems solved in 


i: 
POISDD by the Buneman variant of the CORF (cyclic odd-even 


reduction and factorization) algorithm. Buzbee et al. (1970) 


gtves a complete mathematical description of the algorithm. 
Sweet (1971) gives a brief description of the subroutine. 


A function of two variables, B. j’ is used in the 


numerical model code to represent a two-dimensional array 
Moen provides, on input to POISDD, values of the function 


f. i as well as the specified boundary conditions. On output 
from POISDD, it provides the values of the approximation 
Ë 


T. E For our numerical model, we approximate the outflow 


Boundary condition on input to POTSUD as follows: 


t-at 
* = EN (12) 
ie. ite. ağ 

Po 


On output from POISDD, we correct Eq. (12) by 


EEN 
= + сј 
E? Tı Irax/ C. At 
О 
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